Genome-wide association and Mendelian randomisation analysis among 30,699 Chinese pregnant women identifies novel genetic and molecular risk factors for gestational diabetes and glycaemic traits.

AIMS/HYPOTHESIS: Gestational diabetes mellitus (GDM) is the most common disorder in pregnancy; however, its underlying causes remain obscure. This study aimed to investigate the genetic and molecular risk factors contributing to GDM and glycaemic traits. METHODS: We collected non-invasive prenatal test (NIPT) sequencing data along with four glycaemic and 55 biochemical measurements from 30,699 pregnant women during a 2 year period at Shenzhen Baoan Women's and Children's Hospital in China. Genome-wide association studies (GWAS) were conducted between genotypes derived from NIPTs and GDM diagnosis, baseline glycaemic levels and glycaemic levels after glucose challenges. In total, 3317 women were diagnosed with GDM, while 19,565 served as control participants. The results were replicated using two independent cohorts. Additionally, we performed one-sample Mendelian randomisation to explore potential causal associations between the 55 biochemical measurements and risk of GDM and glycaemic levels. RESULTS: We identified four genetic loci significantly associated with GDM susceptibility. Among these, MTNR1B exhibited the highest significance (rs10830963-G, OR [95% CI] 1.57 [1.45, 1.70], p=4.42×10-29), although its effect on type 2 diabetes was modest. Furthermore, we found 31 genetic loci, including 14 novel loci, that were significantly associated with the four glycaemic traits. The replication rates of these associations with GDM, fasting plasma glucose levels and 0 h, 1 h and 2 h OGTT glucose levels were 4 out of 4, 6 out of 9, 10 out of 11, 5 out of 7 and 4 out of 4, respectively. Mendelian randomisation analysis suggested that a genetically regulated higher lymphocytes percentage and lower white blood cell count, neutrophil percentage and absolute neutrophil count were associated with elevated glucose levels and an increased risk of GDM. CONCLUSIONS/INTERPRETATION: Our findings provide new insights into the genetic basis of GDM and glycaemic traits during pregnancy in an East Asian population and highlight the potential role of inflammatory pathways in the aetiology of GDM and variations in glycaemic levels. DATA AVAILABILITY: Summary statistics for GDM; fasting plasma glucose; 0 h, 1 h and 2h OGTT; and the 55 biomarkers are available in the GWAS Atlas (study accession no.: GVP000001, https://ngdc.cncb.ac.cn/gwas/browse/GVP000001) .

The Born in Guangzhou Cohort Study enables generational genetic discoveries.

Genomic research that targets large-scale, prospective birth cohorts constitutes an essential strategy for understanding the influence of genetics and environment on human health1. Nonetheless, such studies remain scarce, particularly in Asia. Here we present the phase I genome study of the Born in Guangzhou Cohort Study2 (BIGCS), which encompasses the sequencing and analysis of 4,053 Chinese individuals, primarily composed of trios or mother-infant duos residing in South China. Our analysis reveals novel genetic variants, a high-quality reference panel, and fine-scale local genetic structure within BIGCS. Notably, we identify previously unreported East Asian-specific genetic associations with maternal total bile acid, gestational weight gain and infant cord blood traits. Additionally, we observe prevalent age-specific genetic effects on lipid levels in mothers and infants. In an exploratory intergenerational Mendelian randomization analysis, we estimate the maternal putatively causal and fetal genetic effects of seven adult phenotypes on seven fetal growth-related measurements. These findings illuminate the genetic links between maternal and early-life traits in an East Asian population and lay the groundwork for future research into the intricate interplay of genetics, intrauterine exposures and early-life experiences in shaping long-term health.

Verifying the accuracy of self-reported smoking behavior in female volunteer soldiers.

Smoking rates in the military are evaluated through questionnaire surveying. Because the accurate identification of smokers facilitates the provision of smoking cessation services, this study conducted urine cotinine concentration testing to verify the accuracy of self-reported smoking behavior by female volunteer soldiers and analyzed the effects of second-hand smoking on urine cotinine concentrations. This study is a cross-sectional study conducted using purposive sampling on female volunteer soldiers receiving training at the Taichung Recruit Training Center in May 2014. This study simultaneously collected questionnaires and urine samples, and urine samples were analyzed with an enzyme-linked immunosorbent assay. The self-reported smoking rate of female volunteer soldiers was 19.3%, whereas the smoking rate as determined by urine cotinine concentration testing was 26.3%, indicating an overall underestimation of 7.0%. Chi-square (χ2) goodness of fit test results indicated that the distribution of self-reported smoking behaviors and that verified from urine cotinine concentration testing were significantly different. The sensitivity of self-reported smoking behavior was 66.7% with a specificity of 97.6%. There was no significant association between second-hand smoking and urine cotinine concentrations. Questionnaire survey self-reporting methods could underestimate the smoking behavior of female volunteer soldiers and routine testing with biochemical verification is necessary.

CMDB: the comprehensive population genome variation database of China.

A high-quality genome variation database derived from a large-scale population is one of the most important infrastructures for genomics, clinical and translational medicine research. Here, we developed the Chinese Millionome Database (CMDB), a database that contains 9.04 million single nucleotide variants (SNV) with allele frequency information derived from low-coverage (0.06×-0.1×) whole-genome sequencing (WGS) data of 141 431 unrelated healthy Chinese individuals. These individuals were recruited from 31 out of the 34 administrative divisions in China, covering Han and 36 other ethnic minorities. CMDB, housing the WGS data of a multi-ethnic Chinese population featuring wide geographical distribution, has become the most representative and comprehensive Chinese population genome database to date. Researchers can quickly search for variant, gene or genomic regions to obtain the variant information, including mutation basic information, allele frequency, genic annotation and overview of frequencies in global populations. Furthermore, the CMDB also provides information on the association of the variants with a range of phenotypes, including height, BMI, maternal age and twin pregnancy. Based on these data, researchers can conduct meta-analysis of related phenotypes. CMDB is freely available at https://db.cngb.org/cmdb/.

Effect of Surface Morphology Changes on Optical Properties of Silicon Nanowire Arrays.

The optical properties of silicon nanowire arrays (SiNWs) are closely related to surface morphology due to quantum effects and quantum confinement effects of the existing semiconductor nanocrystal. In order to explore the influence of the diameters and distribution density of nanowires on the light absorption in the visible to near infrared band, we report the highly efficient method of multiple replication of versatile homogeneous Au films from porous anodic aluminum oxide (AAO) membranes by ion sputtering as etching catalysts; the monocrystalline silicon is etched along the growth templates in a fixed proportion chemical solution to form homogeneous ordered arrays of different morphology and distributions on the surface. In this system, we demonstrate that the synthesized nanostructure arrays can be tuned to exhibit different optical characteristics in the test wavelength range by adjusting the structural parameters of AAO membranes.

Initial whole-genome sequencing and analysis of the host genetic contribution to COVID-19 severity and susceptibility.

The COVID-19 pandemic has accounted for millions of infections and hundreds of thousand deaths worldwide in a short-time period. The patients demonstrate a great diversity in clinical and laboratory manifestations and disease severity. Nonetheless, little is known about the host genetic contribution to the observed interindividual phenotypic variability. Here, we report the first host genetic study in the Chinese population by deeply sequencing and analyzing 332 COVID-19 patients categorized by varying levels of severity from the Shenzhen Third People's Hospital. Upon a total of 22.2 million genetic variants, we conducted both single-variant and gene-based association tests among five severity groups including asymptomatic, mild, moderate, severe, and critical ill patients after the correction of potential confounding factors. Pedigree analysis suggested a potential monogenic effect of loss of function variants in GOLGA3 and DPP7 for critically ill and asymptomatic disease demonstration. Genome-wide association study suggests the most significant gene locus associated with severity were located in TMEM189-UBE2V1 that involved in the IL-1 signaling pathway. The p.Val197Met missense variant that affects the stability of the TMPRSS2 protein displays a decreasing allele frequency among the severe patients compared to the mild and the general population. We identified that the HLA-A*11:01, B*51:01, and C*14:02 alleles significantly predispose the worst outcome of the patients. This initial genomic study of Chinese patients provides genetic insights into the phenotypic difference among the COVID-19 patient groups and highlighted genes and variants that may help guide targeted efforts in containing the outbreak. Limitations and advantages of the study were also reviewed to guide future international efforts on elucidating the genetic architecture of host-pathogen interaction for COVID-19 and other infectious and complex diseases.

Identifying occult maternal malignancies from 1.93 million pregnant women undergoing noninvasive prenatal screening tests.

PURPOSE: Multiple chromosomal aneuploidies may be associated with maternal malignancies and can cause failure of noninvasive prenatal screening (NIPS) tests. However, multiple chromosomal aneuploidies show poor specificity and selectivity for diagnosing maternal malignancies. METHODS: This multicenter retrospective analysis evaluated 639 pregnant women who tested positive for multiple chromosomal aneuploidies on initial NIPS test between January 2016 and December 2017. Women were assessed using genome profiling of copy-number variations, which was translated to cancer risk using a novel bioinformatics algorithm called the cancer detection pipeline (CDP). Sensitivity, specificity, and positive predictive value (PPV) of diagnosing maternal malignancies were compared for multiple chromosomal aneuploidies, the CDP model, and the combination of CDP and plasma tumor markers. RESULTS: Of the 639 subjects, 41 maternal malignant cancer cases were diagnosed. Multiple chromosomal aneuploidies predicted maternal malignancies with a PPV of 7.6%. Application of the CDP algorithm to women with multiple chromosomal aneuploidies allowed 34 of the 41 (83%) cancer cases to be identified, while excluding 422 of 501 (84.2%) of the false positive cases. Combining the CDP with plasma tumor marker testing gave PPV of 75.0%. CONCLUSION: The CDP algorithm can diagnose occult maternal malignancies with a reasonable PPV in multiple chromosomal aneuploidies-positive pregnant women in NIPS tests. This performance can be further improved by incorporating findings for plasma tumor markers.

Genomic Analyses from Non-invasive Prenatal Testing Reveal Genetic Associations, Patterns of Viral Infections, and Chinese Population History.

We analyze whole-genome sequencing data from 141,431 Chinese women generated for non-invasive prenatal testing (NIPT). We use these data to characterize the population genetic structure and to investigate genetic associations with maternal and infectious traits. We show that the present day distribution of alleles is a function of both ancient migration and very recent population movements. We reveal novel phenotype-genotype associations, including several replicated associations with height and BMI, an association between maternal age and EMB, and between twin pregnancy and NRG1. Finally, we identify a unique pattern of circulating viral DNA in plasma with high prevalence of hepatitis B and other clinically relevant maternal infections. A GWAS for viral infections identifies an exceptionally strong association between integrated herpesvirus 6 and MOV10L1, which affects piwi-interacting RNA (piRNA) processing and PIWI protein function. These findings demonstrate the great value and potential of accumulating NIPT data for worldwide medical and genetic analyses.

Validation of Self-reported Smoking with Urinary Cotinine Levels and Influence of Second-hand Smoke among Conscripts.

Accurate identification of smoking behaviour is crucial to monitor the smoking rate. This study used urinary cotinine (UC) as a biomarker to verify the effectiveness of self-reported smoking behaviour among conscripts during recruit training. The influence of second-hand smoke (SHS) on the UC concentration was also analysed. A cross-sectional study was conducted from July 2014 to December 2014. The participants comprised a total of 621 military service and basic military training conscripts. A self-administered questionnaire survey and a urine test were performed to verify the participants' smoking behaviour. The UC concentration of 100 ng/mL was adopted as the baseline to identify smokers. A high level of consistency was observed between the conscripts' self-reported results and the results validated by the UC concentrations (the overall kappa coefficient was 0.918). Moreover, the overall sensitivity and specificity were 92.9% and 98.1%, respectively. The sensitivity for the military service conscripts was significantly lower than that for the basic military training conscripts (86.1% vs. 97.5%, P-value = 0.002). For the self-reported nonsmokers among the military service conscripts, SHS exposure was related to their UC concentrations. The method of self-reporting through a questionnaire survey can serve as a tool to assess conscripts' smoking behaviour.

Sequencing and de novo assembly of 150 genomes from Denmark as a population reference.

Hundreds of thousands of human genomes are now being sequenced to characterize genetic variation and use this information to augment association mapping studies of complex disorders and other phenotypic traits. Genetic variation is identified mainly by mapping short reads to the reference genome or by performing local assembly. However, these approaches are biased against discovery of structural variants and variation in the more complex parts of the genome. Hence, large-scale de novo assembly is needed. Here we show that it is possible to construct excellent de novo assemblies from high-coverage sequencing with mate-pair libraries extending up to 20 kilobases. We report de novo assemblies of 150 individuals (50 trios) from the GenomeDenmark project. The quality of these assemblies is similar to those obtained using the more expensive long-read technology. We use the assemblies to identify a rich set of structural variants including many novel insertions and demonstrate how this variant catalogue enables further deciphering of known association mapping signals. We leverage the assemblies to provide 100 completely resolved major histocompatibility complex haplotypes and to resolve major parts of the Y chromosome. Our study provides a regional reference genome that we expect will improve the power of future association mapping studies and hence pave the way for precision medicine initiatives, which now are being launched in many countries including Denmark.

Stochastic anomaly of methylome but persistent SRY hypermethylation in disorder of sex development in canine somatic cell nuclear transfer.

Somatic cell nuclear transfer (SCNT) provides an excellent model for studying epigenomic reprogramming during mammalian development. We mapped the whole genome and whole methylome for potential anomalies of mutations or epimutations in SCNT-generated dogs with XY chromosomal sex but complete gonadal dysgenesis, which is classified as 78, XY disorder of sex development (DSD). Whole genome sequencing revealed no potential genomic variations that could explain the pathogenesis of DSD. However, extensive but stochastic anomalies of genome-wide DNA methylation were discovered in these SCNT DSD dogs. Persistent abnormal hypermethylation of the SRY gene was observed together with its down-regulated mRNA and protein expression. Failure of SRY expression due to hypermethylation was further correlated with silencing of a serial of testis determining genes, including SOX9, SF1, SOX8, AMH and DMRT1 in an early embryonic development stage at E34 in the XY(DSD) gonad, and high activation of the female specific genes, including FOXL2, RSPO1, CYP19A1, WNT4, ERα and ERß, after one postnatal year in the ovotestis. Our results demonstrate that incomplete demethylation on the SRY gene is the driving cause of XY(DSD) in these XY DSD dogs, indicating a central role of epigenetic regulation in sex determination.

Discovery, genotyping and characterization of structural variation and novel sequence at single nucleotide resolution from de novo genome assemblies on a population scale.

BACKGROUND: Comprehensive recognition of genomic variation in one individual is important for understanding disease and developing personalized medication and treatment. Many tools based on DNA re-sequencing exist for identification of single nucleotide polymorphisms, small insertions and deletions (indels) as well as large deletions. However, these approaches consistently display a substantial bias against the recovery of complex structural variants and novel sequence in individual genomes and do not provide interpretation information such as the annotation of ancestral state and formation mechanism. FINDINGS: We present a novel approach implemented in a single software package, AsmVar, to discover, genotype and characterize different forms of structural variation and novel sequence from population-scale de novo genome assemblies up to nucleotide resolution. Application of AsmVar to several human de novo genome assemblies captures a wide spectrum of structural variants and novel sequences present in the human population in high sensitivity and specificity. CONCLUSIONS: Our method provides a direct solution for investigating structural variants and novel sequences from de novo genome assemblies, facilitating the construction of population-scale pan-genomes. Our study also highlights the usefulness of the de novo assembly strategy for definition of genome structure.

De novo assembly of a haplotype-resolved human genome.

The human genome is diploid, and knowledge of the variants on each chromosome is important for the interpretation of genomic information. Here we report the assembly of a haplotype-resolved diploid genome without using a reference genome. Our pipeline relies on fosmid pooling together with whole-genome shotgun strategies, based solely on next-generation sequencing and hierarchical assembly methods. We applied our sequencing method to the genome of an Asian individual and generated a 5.15-Gb assembled genome with a haplotype N50 of 484 kb. Our analysis identified previously undetected indels and 7.49 Mb of novel coding sequences that could not be aligned to the human reference genome, which include at least six predicted genes. This haplotype-resolved genome represents the most complete de novo human genome assembly to date. Application of our approach to identify individual haplotype differences should aid in translating genotypes to phenotypes for the development of personalized medicine.

Novel variation and de novo mutation rates in population-wide de novo assembled Danish trios.

Building a population-specific catalogue of single nucleotide variants (SNVs), indels and structural variants (SVs) with frequencies, termed a national pan-genome, is critical for further advancing clinical and public health genetics in large cohorts. Here we report a Danish pan-genome obtained from sequencing 10 trios to high depth (50 × ). We report 536k novel SNVs and 283k novel short indels from mapping approaches and develop a population-wide de novo assembly approach to identify 132k novel indels larger than 10 nucleotides with low false discovery rates. We identify a higher proportion of indels and SVs than previous efforts showing the merits of high coverage and de novo assembly approaches. In addition, we use trio information to identify de novo mutations and use a probabilistic method to provide direct estimates of 1.27e-8 and 1.5e-9 per nucleotide per generation for SNVs and indels, respectively.

Rapid detection of structural variation in a human genome using nanochannel-based genome mapping technology.

BACKGROUND: Structural variants (SVs) are less common than single nucleotide polymorphisms and indels in the population, but collectively account for a significant fraction of genetic polymorphism and diseases. Base pair differences arising from SVs are on a much higher order (>100 fold) than point mutations; however, none of the current detection methods are comprehensive, and currently available methodologies are incapable of providing sufficient resolution and unambiguous information across complex regions in the human genome. To address these challenges, we applied a high-throughput, cost-effective genome mapping technology to comprehensively discover genome-wide SVs and characterize complex regions of the YH genome using long single molecules (>150 kb) in a global fashion. RESULTS: Utilizing nanochannel-based genome mapping technology, we obtained 708 insertions/deletions and 17 inversions larger than 1 kb. Excluding the 59 SVs (54 insertions/deletions, 5 inversions) that overlap with N-base gaps in the reference assembly hg19, 666 non-gap SVs remained, and 396 of them (60%) were verified by paired-end data from whole-genome sequencing-based re-sequencing or de novo assembly sequence from fosmid data. Of the remaining 270 SVs, 260 are insertions and 213 overlap known SVs in the Database of Genomic Variants. Overall, 609 out of 666 (90%) variants were supported by experimental orthogonal methods or historical evidence in public databases. At the same time, genome mapping also provides valuable information for complex regions with haplotypes in a straightforward fashion. In addition, with long single-molecule labeling patterns, exogenous viral sequences were mapped on a whole-genome scale, and sample heterogeneity was analyzed at a new level. CONCLUSION: Our study highlights genome mapping technology as a comprehensive and cost-effective method for detecting structural variation and studying complex regions in the human genome, as well as deciphering viral integration into the host genome.

An atlas of DNA methylomes in porcine adipose and muscle tissues.

It is evident that epigenetic factors, especially DNA methylation, have essential roles in obesity development. Here, using pig as a model, we investigate the systematic association between DNA methylation and obesity. We sample eight variant adipose and two distinct skeletal muscle tissues from three pig breeds living within comparable environments but displaying distinct fat level. We generate 1,381 Gb of sequence data from 180 methylated DNA immunoprecipitation libraries, and provide a genome-wide DNA methylation map as well as a gene expression map for adipose and muscle studies. The analysis shows global similarity and difference among breeds, sexes and anatomic locations, and identifies the differentially methylated regions. The differentially methylated regions in promoters are highly associated with obesity development via expression repression of both known obesity-related genes and novel genes. This comprehensive map provides a solid basis for exploring epigenetic mechanisms of adipose deposition and muscle growth.

Bisulfite sequencing reveals that Aspergillus flavus holds a hollow in DNA methylation.

Aspergillus flavus first gained scientific attention for its production of aflatoxin. The underlying regulation of aflatoxin biosynthesis has been serving as a theoretical model for biosynthesis of other microbial secondary metabolites. Nevertheless, for several decades, the DNA methylation status, one of the important epigenomic modifications involved in gene regulation, in A. flavus remains to be controversial. Here, we applied bisulfite sequencing in conjunction with a biological replicate strategy to investigate the DNA methylation profiling of A. flavus genome. Both the bisulfite sequencing data and the methylome comparisons with other fungi confirm that the DNA methylation level of this fungus is negligible. Further investigation into the DNA methyltransferase of Aspergillus uncovers its close relationship with RID-like enzymes as well as its divergence with the methyltransferase of species with validated DNA methylation. The lack of repeat contents of the A. flavus' genome and the high RIP-index of the small amount of remanent repeat potentially support our speculation that DNA methylation may be absent in A. flavus or that it may possess de novo DNA methylation which occurs very transiently during the obscure sexual stage of this fungal species. This work contributes to our understanding on the DNA methylation status of A. flavus, as well as reinforces our views on the DNA methylation in fungal species. In addition, our strategy of applying bisulfite sequencing to DNA methylation detection in species with low DNA methylation may serve as a reference for later scientific investigations in other hypomethylated species.

Systematic assessment of reduced representation bisulfite sequencing to human blood samples: A promising method for large-sample-scale epigenomic studies.

Complementary to the time- and cost-intensive direct bisulfite sequencing, we applied reduced representation bisulfite sequencing (RRBS) to the human peripheral blood mononuclear cells (PBMC) from YH, the Asian individual whose genome and epigenome has been deciphered in the YH project and systematically assessed the genomic coverage, coverage depth and reproducibility of this technology as well as the concordance of DNA methylation levels measured by RRBS and direct bisulfite sequencing for the detected CpG sites. Our result suggests that RRBS can cover more than half of CpG islands and promoter regions with a good coverage depth and the proportion of the CpG sites covered by the biological replicates reaches 80-90%, indicating good reproducibility. Given a smaller data quantity, RRBS enjoys much better coverage depth than direct bisulfite sequencing and the concordance of DNA methylation levels between the two methods is high. It can be concluded that RRBS is a time and cost-effective sequencing method for unbiased DNA methylation profiling of CpG islands and promoter regions in a genome-wide scale and it is the method of choice to assay certain genomic regions for multiple samples in a rapid way.

Structural variation in two human genomes mapped at single-nucleotide resolution by whole genome de novo assembly.

Here we use whole-genome de novo assembly of second-generation sequencing reads to map structural variation (SV) in an Asian genome and an African genome. Our approach identifies small- and intermediate-size homozygous variants (1-50 kb) including insertions, deletions, inversions and their precise breakpoints, and in contrast to other methods, can resolve complex rearrangements. In total, we identified 277,243 SVs ranging in length from 1-23 kb. Validation using computational and experimental methods suggests that we achieve overall <6% false-positive rate and <10% false-negative rate in genomic regions that can be assembled, which outperforms other methods. Analysis of the SVs in the genomes of 106 individuals sequenced as part of the 1000 Genomes Project suggests that SVs account for a greater fraction of the diversity between individuals than do single-nucleotide polymorphisms (SNPs). These findings demonstrate that whole-genome de novo assembly is a feasible approach to deriving more comprehensive maps of genetic variation.